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Abstract 

We give an upper bound in 0(d' n+1 ^ 2 ) for the number of critical 
points of a normal random polynomial with degree at most d and n vari- 
ables. Using the large deviation principle for the spectral value of large 
random matrices we obtain the bound 



O (exp(-/?n 2 + ~ log(d - 1))) 



(/3 is a positive constant independent on n and d) for the number of 
minima of such a polynomial. This proves that most normal random 
polynomials of fixed degree have only saddle points. Finally, we give a 
closed form expression for the number of maxima (resp. minima) of a 
random univariate polynomial, in terms of hypergeometric functions. 

1 Introduction 

We consider a random polynomial / over the reals with n > 1 variables and 
degree at most d > 2. The problem is to compute, on the average, its number 
of critical points (the number of real roots of the system Df(x) — 0), and its 
number of local minima. Since a generic polynomial has only nondcgcncratc 
stationnary points, this last number is also given by the real roots of the system 
Df(x) = such that D 2 f(x) is positive definite. This reduces our problem 
to the computation of the number of real roots of a polynomial system under 
certain constraints. 
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Generally speaking, let F = (Fi, . . . ,F n ) be a random system of real poly- 
nomial equations with n variables and degree Fj < dj. Let N F (U) denotes 
the number of zeros of the system F{x) = lying in the subset U C R™ and 
N F {W L ) = N F . Little is known on the distribution of the random variable 
N F (U) . A classical result in the case of one polynomial of one variable is given 
by Kac [3] , [TU] , who gives the asymptotic value 

E(N F ) « -logd 

as d tends to infinity when the coefficients of F are Gaussian centered indepen- 
dent random variables with variances equal to 1. But, when the variance of the 
i— th coefficient is equal to (.) (Weyl's distribution), we have (see Bogomolny- 
Bohias-Leboeuf [5] and also Edelman-Kostlan [6]) 

E(N F ) = Vd. 

In 1992, Shub and Smale extended this result to a real polynomial system 
F where 

Fi{x\, . . . , X n ) ^ ^ (Xi^ a X^ . . . X n , 

ai + . . .+a n <_di 

when the coefficients ai jQ are Gaussian centered independent random variables 
with variances equal to 

(d\ = ck\ 

(see Kostlan [11] on this distribution and its properties). Their result is 

E(N F ) = y/di... d n 

that is the square root of the Bezout number of the system. 

A general formula for the expected value of N F (U) when the random func- 
tions Fi, 1 < i < n, are stochastly independent and their law is centered and 
invariant under the isometries of R™ can be found in Azai's-Wschebor [3] . This 
includes the Shub-Smale formula as a special case. 

This result has also been extended by Rojas [H] to multi-homogeneous poly- 
nomial systems, and then partially by Malajovich and Rojas [12] to sparse poly- 
nomial systems. 

Wschcbor in [T7] studies the moments of N F and Armentano-Wschcbor [5] 
consider random systems of equations of the type Pi(x) + Xi(x), 1 < i < n, 
x € R n , where the P/s are non-random polynomials (the signal) and the X[s 
are independent Gaussian random variables (the noise). 

Notice a major difference between these studies and the case considered here: 
the n equations of the system Df{x) = are not independent! 

Through this paper we denote by V = Vd.n the space of degree at most 
d, n-variate polynomials with real coefficients. This space is endowed with the 
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inner product: 

if'9)v = J2 [a] fa9 ° 

\a\<d ^ ' 

where a = (aci, . . . , a n ) € N™ is a multi-integer, \a\ = a\ + 

|a|<d |a|<d 

and again 

'(A d! 



ai!...a„!(d- |a|)! 
We make a probability space in considering the probability measure 

1/2 



i.e. a random polynomial has here Gaussian centered independent random co- 
efficients with variances equal to 

Let S n be the space ofnxn real symmetric matrices, endowed with the 
Frobenius inner product (R, S) = Trace(i? T S') and its induced norm 

\\sr= £ sf r 

l<i ;< 7<n 

The Gaussian Orthogonal Ensemble is the space S n together with the probability 
measure 

e -||S|| 2 /2 e -\\S\\ 2 /2 



)/i dS 2"/ 2 7r"(™+ 1 )/ 4 A dSi3 ' 
l<i<j<n 

Thus, the entries of a matrix in iS„ are independent Gaussian random variables 
with mean and variance 1 for a diagonal entry, and mean and variance 1/2 
for a non-diagonal entry. 

Our first main result is the following: 

Theorem 1. Let Cd, n denote the expected number of critical points of a random 
polynomial of degree at most d in n variables, and E dn the expected number of 
minima. Let P n be the probability that a matrix in the Gaussian Orthogonal 
Ensemble is positive definite. Then, for every n > 2, 

C2, n = 1 and E 2 ,n = Pn, 

and for d > 3 

C d , n < J^(d - l)(" +2 >/ 2 and E d . n < J^(d l)("+ 2 )/ 2 P n . 
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When n = 1 one has 



C dj i = 2S d ,i = / ,, 2± .w, ±1 , air < 1 + - 2. 

7r 7 (c«H + l)(r- + 1) 



Moreover, when d — > oo, 

C<i,l 



1. 



1 + 2 

Let P n be the probability that a matrix in the Gaussian Orthogonal Ensem- 
ble GOE{n) is positive definite: 

P " ~ /-++ 2»/2 7r n(»+l)/4 A 

" l<*<j<n 

Via the change of variable 5 = QAQ T with Q S O n and A = diag(Ai > . . . > 
A n > 0) one has 



2" Jnii<- (27r) n (" +1 )/ 4 



Vol 



where A € R™ if and only if Ai > . . . > A n > and 

2 n(n+3)/4 r Q / 2 )n(n+l)/2 

= n;=ir((n-j + i)/2) 

(see Mehta [T3] for the description of P n as an integral over M. n and Federer [7] 
for the volume of the orthogonal group) . The following values are easy to obtain 

1 2-V2 7T-2V2 

Pi = -, Pi = — , P3 = — • 

2' 1 4 ' J 4tt 

P3 was computed by Carlos Beltran. 

Using the large deviation principle for the spectral value of large random 
matrices (see the appendix at the end of this paper) we see that the asymptotic 
value of P n for large values of n satisfies 

limsup log(P„) < -a 

n — y 00 ri 

where a is a positive constant independent on 11. Thus, there exist two positive 
constants (3 and 7 such that, for every n > 1, 

Pn <ie- 0n2 . 
This gives our second main theorem: 

Theorem 2. There exist two positive constants (3 and K such that for every n 
and d the number of minima of a random polynomial satisfies 

E d . n < Ke-> 3n2+2 i los( - d - 1) . 
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Remark 1. This is a quite surprising result : it shows that most of random 
polynomials of reasonable degree have only saddle points. 

We want to thank here Alice Guionnet, Manjunath Krishnapur and Balint 
Virag who introduced us in the world of large deviation delicacies. 

2 The space of n-variate polynomials 

The inner product space V,{-,-)-p has several interesting properties resumed in 
the following 

Lemma 1. 1. It admits the reproducing kernel K(z,x) = (1 + (z,x)) d : 

f(x) = (K(.,x)J) v (1) 
for any x G W l and f G V . 

2. It has a representation formula for the derivatives: for any integer k > 1 
and x, Ui, ■ ■ ■ ,u k € M™ we have 

D k f(x)(u!,--- ,Uk) = (K k (.,x,ui,--- ,u k ),f) v , (2) 

with 

K k (z, x, mi, ■ • • , tifc) = D k K(z, x)(ui,--- , u k ) = (3) 
d ■ ■ ■ (d, — k + 1) (z, U\) ■ ■ ■ (z, u k ) (1 + (z, x)) d k . 

3. This scalar product is orthogonally invariant: 

(foU,goU) v = (f,g) v (4) 
for any f,g(zV and the orthogonal transformation U G O n . 

Proof. The first two formulas are well known and easily obtained via a direct 
computation. For the orthogonal invariance see [1], section 12.1, or [TT]. □ 

A second interest of Wcyl's distribution for polynomials is due to the fol- 
lowing identity: let f(x) = x T Sx (here S is a symmetric n x n matrix) be a 
homogeneous degree 2 polynomial, then ||/||-p = 11511- This is the reason why 

Proposition 1. C2, n = 1 and E2, n = Pn- 

Proof. Since a generic degree 2 polynomial has only one critical point we have 
C2, n = 1- Given / G V2, n we can write it 

f(x)=a+ 2J h%i+ ^2 a a x i + a^XiXj. 

l<i<n l<i<n l<.i<j<n 

One has 

ii/ii^ = « 2 + ^ E b "+ E ^ + \ E 4 
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so that 

J D'f (o)>o 2™(™+ 1 )/4(2^)(«+ 1 )("+ 2 )/ 4 



D- 



e -(E I a?, + iE I<3 a%)/2 

m>0 2"(«-i)/4(2^)«(«+i)/4 da - 
i n2 , 



To compute this last integral we let S = i£> 2 /(0); this gives 

r P -lls|| 2 /2 

rp . / nn(n-l)/2 jc _ p 

~ /s>o 2«(«- 1 )/ 4 (2^)™("+ 1 )/ 4 



□ 



3 An integral formulation 

Let us define 

evali : V x E" -> R", cvali(/,x) = D/(x). 
The incidence variety for real critical points of a polynomial is defined by 

V = {(f,x) ePxR" : evali(i» = 0} . 

The derivative of evah is given by 

Devali(/,a;)(/,i) = Df(x) + D 2 f{x)x 

for any /,/£? and £, i; G R™. Since this derivative is onto, is a submanifold 
and its dimension is 

dim V = dim V = ( 

The tangent space at (/, x) £ V is given by 

T (/)X) y = ker J Devali(/ ) x) = {(/,*) eV xl» : Dj(x) + D 2 f(x)x = j . 

The restriction tt2 : V — > IK" of the projection P x R" — > R" is surjective 
and is also a regular map because for any (/, x) <E V the derivative DiT2(f,x) : 
T( ftX) V -» R™ is surjective. The fiber of tt 2 above x € R" 

V x = {(/, x) G P x R" : evali(/, x) = 0} 

is isomorphic to a dim V — n linear space. V x is equipped with the volume form 
inherited from the induced metric. 

The restriction tti : V — ► "P of the projection "P x R" — > P is a smooth map. 
A given / G "P is a regular value of ~k\ when either / has no critical point or 
when, for any x such that (/, x) 6 V, Dit\(f, x) : T^^V — > V is surjective. This 
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last condition is satisfied when the second derivative D 2 f(x) is an isomorphism 
which is the generic situation: 

E' = {(f,x)€V : dct D 2 f(x) =0} 

is a submanifold in V and dimS' < dimV^. Thus £' and its image £ = 7Ti(S') 
have zero measure and we may ignore them. For any (/, x) € V \ £' and any 
/ e V we have Dm(f,x)(f,x) = f for x = -D 2 f{x)- 1 Df{x) and the fiber 
above / 

V f = {(f,x)eVxR n : evah(f,x)=0} 

consists in a finite number of points. 

Given (/, x) € V \ £' we are in the context of the implicit function theorem 
that is V is locally around (/, x) the graph of the function 

G = 7T2 O TTi 1 

where n^ 1 is the local inverse of tt\ such that 7r^ 1 (/) = (f,x). Since the graph 
of DG(f) is the tangent space JVy _ X )V we get 

DG(/)/ = -LPfix^Dfo) (5) 

for any / e ?. 

Like in [4] section 13.2, theorem 3, we have the following 

Proposition 2. Let U be a measurable subset of V . Let us denote by #(/, U) 
the number of pairs (f,x) G U and by Ejj the expectation of #(f,U) when f is 
taken at random: 

With these notations, one has 

r r , e-w^i 2 

E V dx dct(DG(f)DG(fyr 1/2 (0 , dimV/2 dV x . (7) 

Remark 2. In our context two sets are of particular interest: U = V to compute 
the average number of critical points of a polynomial Cd, n , and U = V+ with 

V+ = {(f,x) EV xR n : Df(x) = and D 2 f(x) > 0} 

(here > means positive definite) for the average number of local minima Ed, n - 

We have now to compute the determinant appearing in equation [7] This is 
done in the following 

Proposition 3. Under the notations above 
dct{DG(f)DG(f)*) = d n {l + Wxf)^- 1 ^- 1 ^ + d\\x\\ 2 ) \detD 2 f(x)\~ 2 . (8) 
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Proof. Let us denote Df(x) = D x f. Since DG(f)f = -D 2 f(x)- 1 D x f and 
since D 2 f(x) is symmetric, we get 

DG(f)DG(fy = D 2 f{x)- 1 D x D%D 2 f{x)- 1 

so that 

det(£>G(/)£>G(/)*) = dct(D x D*) |det D 2 f(x)\~ 2 . (9) 

To compute det(D x D*) we use the representation formula for the derivative 
(equation [2]) with k = 1. Let us denote by ej, 1 < i < n, the canonical basis in 
E". Then, for any /eP, 



i 

so that, with x S R™, x = Xjej, 
D* x ij) = (x,D x f\ (x,^2e i /K 1 (.,x,e i ),f\ ) ^ x, x, e 4 ), / 



Thus, we get 



= '^2iiiKi{.,x,e i ) 

i 

and consequently 

D x D* x x = y^ej( K 1 (.,x,e i ) 7 ^2x j K 1 (.,x,e j ) 



Xj-a [z, ej) (L ■+■ [z, x) ) 

3 



d{d-l)x i x j {l + \\x\\ 2 ) d -' 2 iH^j 



d(d - l)x 2 (i + ||x|| 2 ) d - 2 + d(l + Wxfy- 1 if i = j 

" 'J 

which correspond to the matrix 

d(d - 1)(1 + ||x|| 2 ) d ~ 2 xx T + d(l + \\xf) d - x I n . 

Its eigenvectors are x and any nonzero vector in the orthogonal subspace x^. 
The corresponding eigenvalues are 

d(d 1)(1 + ||x|| 2 ) d - 2 ||x|| 2 + d(l + Hxll 2 )^ 1 = d(l + ||x!| 2 ) d - 2 (l + d||x|| 2 ) 

with multiplicity 1, and 

dCl + Hxll 2 )^ 1 

with multiplicity n — 1 so that 

det D X D* = d"(l + llxll 2 )™^- 1 )- 1 ^ + d||x|| 2 ). 
Our proposition combines this value and equation [51 □ 



If we combine propositions [2] and [3] we obtain the following integral formu- 
lation 

Proposition 4. Let U be a measurable subset ofV. One has 

E = f f \detD 2 f(x)\ er\\nl/2 

(10) 

An action of the orthogonal group O n on V x M" is defined by 

(Q, /, x) £ O n xPxl"^ (/ o Q, Q T x) ePx M™. 

This action leaves the incidence variety V invariant and also the scalar product 
(., .) v (lcmma[T]). For this reason, when the measurable set U is itself invariant, 
the integral on V x n U in proposition 2] only depends on r = \\x\\. Thus, taking 
spherical coordinates in M. n , we get: 

Proposition 5. Let U be a measurable subset ofV invariant under the action 
of O n (for any (Q,f,x) £ O n x U we have (/ o Q,Q T x) £ U). Under this 
condition 



q„ f°° r n ^dr f 2 e -\\f\M* 

d^Jo RV-^-^dri + iy/ 2 J Vre nu ^ (27r) dim W2 ' 



where a„ = ^°^„ /2 = 2 „ /2r 2 (ra/2) , -R = vV 2 + 1 and ref = (r, 0, . . . , 0) . 

Remark 3. The measurable sets considered here: [7 = 1/ and [7 = V" + = 
{(/, x) £ V : D 2 f(x) > 0}, are clearly invariant under the action of O ra . 

4 The inner integral 

Our objective is now to compute the integral over V rei n U appearing in propo- 
sition o 

Let D 2 : V rei — > S n denote the operator / i— > D 2 f(rei). We would like 
to compute its pseudo-inverse : S n — > (kcrZ? 2 )- 1 -. This means that is the 
minimum norm right inverse of D 2 (D 2 = ids„ ) . 

This will allow us to 'integrate out' ker/J 2 : 

-Il/ll 2 / 2 

|det£> 2 /| 6 ^. T ^ /0 dV; ei = (11) 

/ IdetgHdet^^l 1 / 2 2 d5. 

JD 2 (v rei nU) ( Zn ) ' 

To compute ^(S) and | det \P* \I/ 1 we need the following lemma: 

Lemma 2. Let us denote 



• ej, 1 < i < n, the canonical basis in R™, 

• d Ct = K 1 (z,re 1 ,e l ), 

• d eze] = K 2 (z,rei,e i ,e 1 ), 

• R = s/l + r 2 . 
Then, 

1. (d ei ,d ei ) p d(l + dr 2 )R 2d - i 

2. Ifi^l, then (d ei ,d ei ) v = dR 2d - 2 

3. Ifi^j, then(d ei ,d ej ) v =0 

I (d ei ,d eiei ) v ^d{d^l)(dr 2 + 2)rR 2d - (i 

5. If(i,j,k) ? (1,1,1), then {d e] ,d etek ) v =0 

6. (d eiei , d eiei ) v = d{d - 1) (d(d - l)r 4 + 4(d - l)r 2 + 2) i? 2d - 8 

7. Ifk^l, then (d eiek ,d eiek ) v =d(d-l){{d-iy + l)R 2d - 6 

8. Ifi £ 1 and fc ^ 1 , then {d eiek , d^ h ) v = (1 + fo) d(d - l)i? 2d - 4 (S ik is 
the Kronecker symbol), 

9. If{i,k} ^ {j,l}, then (d eiCk , d ejei ) p = 

Proof. It is a consequence of the representation formulas given in lemma [T] 

• (d ei ,d ei ) v = (Ki(.,re 1 ,ei),I<i(.,rei,e 1 )) v = ^■Ki(z,re 1 ,ei)\ z = rei = 
£dzi(l + rz^ 1 | z=rei = d(l + r 2 ) d ~ 2 (l + dr 2 ), 

and similarly 

• {9ei,d ei } v = ■£ T Ki(z,re 1 ,e i ) | a=rei = ^-dz,{l + rz^' 1 \ z=rei = d(l + 
r 2 ) d ~\ 

• (de i: d ej ) v = ^K 1 (z,re 1 ,ej)\ t = rei = -§^dz (\ + rzi) d_1 | 2=rei = 
when i j, 

• {d ei ,d eiei ) T> = ^■K 2 {z,re 1 ,e 1 , e t ) \ z = rei = -^d(d-l)zf(l+rzi) d ~ 2 | 2=rei = 
d(d- l)r(2 + dr 2 )(l + r 2 ) d - 3 , 

• {de ] -,d eiek ) v = -£jd(d- l)ziZ k (l + rz x ) d ~ 2 \ z=rei = when (i,j,fc) 7^ 
(1,1,1), 

. (d eiei ,d eiei ) v = -^d(d-l)z 2 (l + r Zl ) d - 2 \ z=rei = d(d-l)(l + r 2 ) d - 4 (2 + 
4(d-l)r 2 + d(d-l)r 4 ), 
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(d eiek ,d eiek ) v = ^—d(d-l)z 1 z k (l + rz 1 ) d - 2 \^ rei = d(d-l)(l + 
r 2 ) d - 3 (l + {d-l)r 2 ), 



(d ei e k ,d ei e k ) v = -£—d(d-l)z l z k (l + rzi) d 2 \ z=rei = (1 + S ik )d(d 
\){l + r 2 ) d - 2 , 



{d ei e k ,d ejei ) v = gr-^d(d- l)zjZi(l + rzi) d 2 \ z=rei = when {i,k} ^ 



□ 

Let us now evaluate 4'. Recall that 

V rei ={f€V : D/(rei) = 0} 

or, in other words, / € V rei if and only if 

{f,d ei ) r =0, l<i<n. 

Thus, by lemma [21-3, d ei , 1 < i < n, constitue an orthogonal basis of . We 
also have 

ker D 2 = Span {d e%ej , 1 < % < j < n} X H V rei 

hence, 

(kerZ? 2 )- 1 = Span {Pd eiej , 1 < i < j < n} 

where P stands for the orthogonal projection onto V rei . We have seen that for 
{i,j,k) ^ (1, 1, 1), deiCj _L de k (lcmma[2]-5). Hence, 

£ 161 6 161 u ei mo n2 

\\oei\\ v 

and for ^ (1,1), 

Pd, 



Let us now show that 

Pd„ 



*(S)= X! s 



l<i<j<n 11 e > e j"P 



Since this expression is clearly in (ker D 2 )- 1 it suffices to prove that D 2 o ^(S) 
S for any S € 6>„ i.e. 

DH{S){re l ){e k ,e l )=S u 

or, using lemma [1} that 



(de k e n Sy V 3 / ~~ Ski- 

\ 1<«<J<™ II ^ e i e i\\v I -p 

This last equality holds because Pd eiej , 1 < i < j < n, constitue an orthogonal 
basis of (kerZ? 2 )- 1 . 
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It is important to have in mind that ^ is not an isometry, we have 



l<i<j<n \\ r0 ^i\\V 



We introduce now the functions 
A(d,r)-- 



d(d - l)r 4 + 2dr 2 + 2 
{dr 2 + l)i? 4 



and 



B(d,r) = 



(d- l)r 2 + 1 



R 2 



where again R = y/1 + r 2 . 
Lemma 3. Let i < j. Then, 

{A(d, r) 2 ifi = l and j = 1 

B{d,r) 2 ifi = l andj ^ 1 

(1 + 5ij) if i 1 and j 7^ 1 

with Sij = 1 when i = j and otherwise. 

Let us now compute det For any / = Yli<i<j< n fn Pd < ■ e (ker D 2 )- 

and for any S € S n we have 



(**(/),S) = (/,*(5)), 



Therefore, we have always for any T £ S n : 



!<2<j<n 



<***(T),S) = £ 



J 1 , j 



\<i<j<n 1 1 P9 ei e . J J 



We write the matrix of the operator \&* ^ with respect to the orthonormal basis 
of S given by eief , . . . ,e n e^ and then, for i < j, ^= (e,ej + e 3 ef )): 



2||pa ci e 2 p 



Using lemma Owe obtain: 
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Lemma 4. 

(det vr *)V2 = ( d ( d _ l)i? 2d - 4 )"^ A(d, r)- 1 ^^ r)^"- 1 ) . 

At this point 
Proposition 6. Under the conditions above, 



d n/2 J q ( dr 2 + 1) l/2 i? (d-l)„-l J D2{UnVriii) ( 27r )dimS„/2 

In particular, 



On_ l~ (det VP**)* r^dr /■ JetSL_ |t(s) |y 2 ,„ 
rf "d«/ 2 7 (dr 2 + l)l/2i?(d-Dn-l (2 ^)dim5„/2 e 



and 



'd™/ 2 J (dr 2 + l)V2i?(d-l)"-l 7^++ ( 27r )dim5„/2 

where S„ + denotes the set of positive definite matrices. When n = 1 , 



Q,l - 2^,1 - — — ^ {dr 2 + 1){r 2 + 1) d '- 

Proof. The three first formulas are obtained in combining proposition equa- 
tion [Tl] and lemma 31 For the case n = 1 we obtain 

2 Z" 00 dr f°° s s 2 

f], , — / / p 2d(d-l)H 2d - 4 A 2 r| o = 

dA dVd~JV2^Jo A(dr 2 + l)V2ijM-4 7 V2T 



Vd - 1 Z" 00 */d(d - l)r 4 + 2dr 2 + 2 , 

-dr. 



(dr 2 + l)(r 2 + 1) 

The identity Cd,i = 2.Ed,i is easy. □ 

5 Some integral lemmas 

The term e II II / 2 in the inner integrals of Proposition [5] can be simplified 
through additional changes of coordinates. We reparametrize the spaces S n and 
S++ though a stretching S^T = A^SA" 1 . 

The stretching coefficients arc A; = (2d(d- l)i? 2d ~ 4 ) 1/4 for i > 2, A : = 
B(d,r)A 2 and, A = Diag(Ai, A 2 , . . . ,A„). We obtain 



_}_ [ $11 , \ " , \ " 1 c-2 ] 

l)i? 2d - 4 \ A 2 ^ B 2 ^ 1 + Su ij 

1 \ j=2 l<i<j<n 3 ) 
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and 
so that 



J f J^n + V ^ + V 1 s* 2 I 

l)i? 2d - 4 I 2B A B 2 ^ 1 + <fo ij 

' \ j=2 Ki<j<n lJ J 



||*(S)||?, = ||a- 1 sa- 1 || 2 + 

Let us define T = A^SA' 1 . We get 



1 1 



S - 



2 



A 2 2B A J d{d-l)R 2d - 



so that, via this change of variable, 



d 2 (U) V2tt 



\detS\ e -||*(s)|| 2 /2 dS - 



n+3 



--1 IT, 



nAj;) I M e - ir nw T «;dr. 



If C/ C V, we define the auxiliary quantity 

Cu{d,r,n)= I ^ dimSn e V V 7 ;dT - 



A- 1 D 2 (unV r )A- 1 y/2ir 

There are two cases of interest corresponding to U = V for the average of 
critical points and U = V+ for the average number of local minima. The corre- 
sponding functions are denoted Cy(d, r, n) and CV + (d, r, n). Using proposition 
[5] we get (the proof is easy and left to the reader) 

Proposition 7. 

_ 2V2(d-l) n / 2 [°° {(d-l)r 2 + l) 2 r"- 1 

E v — — - — / / ===== -Cu(d,r,n)dr 

r(n/2) J R 2 ^d(d-l)r± + 2dr 2 + 2R n - 1 V ' 

Moreover 

JdetTl 

S„ V27T d 



CV(d,r,») = / _ dim ^ e ^V" J -J AT 



and 
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6 Proof of Theorem [T] 

To prove our main theorem we use both the proposition [7] and the case d = 2 
already investigated in the proposition [1] We have 

1 n 2V2 [°° r n ~ x C v (2,r,n) 

1 = C 2 ,n = —, — I -7T—T 7= dr 



T(n/2)J V2 
and 

Cv%r,n)= I , -iflmiWi*)^ 
Js n (27r) ,l (" +1 )/ 4 

Lemma 5. The quantity A(d, r) = 2 ^^'^ 2 — 1 satisfies, for all r > and d > 2, 
£/ie scaling law: 

X R , 

A(2, rVd~" 1) < A(d, r) < A(2, -yrA/d - !) 



Proof. We write 



. . . , , . , d-2 (d- l)r 
A(d,r) = 2(d- l)r 2 + 



d (d - l)r 4 + 2r 2 + § ' 
The lower bound is now obvious. The upper bound is obtained as follows: 

»/ , x , n 2 - 2 (d- l)r 4 

A(d,r) = 2(d- l)r 2 1 



d (d- l)r 4 + 2r 2 + 



ii 



< A(d,r) = 2(d- l)r 2 + ^-^(d- l)r 2 

< ^A(2,r 2 Vd~ T) 
V5 



A(2,^r 2 Vd~T). 



□ 



It follows from Lemma [5] that 



C v (d,r,n) = / J^e-Kimi^)^) 



/— — dim £"„ 
S„ V27T 



< 



|detT| c _i(|[Tii 2 +A(2,rVd^Dr 1 2 ,) 

1— — dim Sn 

C v (2,rVd-l, n) 



and similarly Cy+(d,r,n) < Cy + (2, r\Jd — 1, n). Now we have: 

2V2(d-l)"/ 2 /•<» ((d-l)r 2 + l) 2 r- 1 

Gj„ = ; — — / T Cv(d,r,n) dr 

r(n/2) 7o i? 2 v /d(d-l)r 4 + 2dr 2 + 2i?"- 1 V ' 



2V2(d-l)"/ 2 f 00 ((d-l)r 2 + l) 



2 



r(n/2) 7 i?2^rf(d-l)r 4 + 2dr 2 + 2i?"- 1 V ' 
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We set s = rV d — 1 and S = \/d— 1 + s 2 to obtain: 

2^(d-i)»/ 2 r (rf-i)^ 2 _ , 9 . i 



r(n/2) y 5 2^ S 4 + 2 _4_ S 2 + 2 5«-i ' v^i 

Since 

(d-i), 2 < . /rf - T 



and C2, n = 1 we obtain 

(d-l)(»+2)/2^/2 

Gd > n ~ Td 

and the same argument holds for Ed n . 



7 The Riemann surface 

We rewrite the case n = 1 (proposition [6|) for convenience as: 

E d , = / 9 (z)dz (12) 



2tt 



with 

sO) = 



^ 4 + ^ + 



{I + z 2 ){\ + dz 2 ) 

At this point we encounter a classical situation: we want to compute a 
line integral of a function g(z), which is a two- branched meromorphic function 
of C. In order to apply the residue theorem, we need first to replace g by a 
regular meromorphic function, defined in the relevant Riemann surface R. The 
branching points of the Riemann surface are the roots of the polynomial inside 
the square root. If we set 



< 



-1 + ul 



d-l 



with the branch of the external square root in such a way that £ belongs to the 
positive quadrant, we can now factorize 

, 4 + ^ 2 + ^(W)(,-C> + C)(, + C). 

It follows that the Riemann Surface R is a twofold cover of C with branch 
points C, -C, -(, (■ 
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Figure 1: The Linear Fractional Map w z = ty(w) 



Let 7 be the arc of circle (centered in the origin) joining — £ to £ crossing 
the positive imaginary axis. Notice that it croses the segment [i / y/d, i] . Let T> 
denote the upper half plane with 7 removed. 

Then, the positive branch of ^/z 4 + j^z 2 + 2{d-i) on ^ extends to a unique 
branch on T>. The square root is real and positive on [0, £|£|] and real and 
negative on [i|£|,ioo). 

The residue theorem is now: 



g(z)dz-2 / g(z)dz2iriRes [z=i/V j ] g(z) + 2TriRcs [z=i] g(z) 

Residues are respectively , ~\ _ and ~ i y d ~ 2 - . Therefore, 

1 Vrf - ^ (d-l)Vd f , , , 
^,i=2+^— + - ^J— j g{z)dz (13) 

(We mean the integral of the branch that is positive on i|£|). 
Now, in order to integrate g(z), we introduce a linear fractional transforma- 
tion mapping the real line onto the circle containing 7. Namely, 

T , , Aw + B 

* H = c^Td 
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with A = |C|, B = i\C\, C = i,D = l. For the record, AD - BC = 2|C| 



Let s 



_ Rc(C) 



Define also si 



i-KI „ 
T+ICT' 32 



KI+Im(C)' 

S4 = S3 . We have the following mapping table for <5: 



l+ICI V3 



and 









^(w) 


-1 


-ICI 


is\ 


i 





i|CI 


is 2 


—i 


1 


ICI 


is 3 


i/Vd 


-s- 1 


-c 


iS4 


-i/Vd 


— s 


-c 






s 


c 






s" 1 


c 







Changing coordinates, 

g(z)dz = 2c(d)Re 



y/ (w 2 — S 2 )(w 2 



dw 



with 



c(d) 



(AD - BC)JA 4 



^A 2 C 2 



2 r 4 



(A 2 + C 2 )(dA 2 + C 2 ) 



G 0(cT 3/2 ) 



(More precisely: hW 3 / 2 C«) = -2 7 / 4 V|-^ ~ -6.2151). 
At this point, good practice seems to be: 

1. Multiply numerator and denominator by the conjugate of the denomina- 
tor, in order to obtain a real polynomial in the denominator. 

2. Multiply numerator and denominator by the square root. 

3. Expand in partial fractions. 

4. Put into Legendre normal form. 

5. Write down the integral in terms of elliptic functions K and II. 
We expand the integrand in partial fractions: 



g(z)dz 



2c(d) 
2c(d) 
2c(d) 



1 



vV 2 


— s 2 )(w 2 — s 2 ) 




1 




— s 2 )(w 2 — s~ 2 ) 




1 



[0,s] y/(w 2 - S 2 )(w 2 - S 2 ) 



l + J]Re 



Rk 



fe=i 



W - ISA; 



dw 



1 I V s k 2 Rc(R k (w + is k )) 



1 + 



k -l + w 2 s k \ 



dw 
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Pole 


Residue R k 


Argument rife 


isi 




(s 1 -s 2 )(s 1 -s 2 




_ s 2 


i 


0l-s 2 )(si-s 3 )0i- 


S 4 ) 




is 2 




(S2-S 2 )(S2S- 2 






i 


{S2-Sl){s 2 -S 3 )(s 2 - 


S4) 




is 3 




(s 3 s 2 )(s 3 - S - 2 , 




2 

S 


i 


(S3-Sl)(«3-S2)(s 3 - 


s 4 ) 


_ P" 

s 3 


is4 




(s 4 -s 2 )(s 4 -s~ 2 




2 

S 


i 


(S4— Sl)(s4— S2)(S4~ 


«3) 





Table 1: Residues and arguments 



(the last step uses the fact that all residues R k are pure imaginary). Residues 
are given in Table [1] We use formula [1] [17.4.45] to compute the parameter 
m = ,s 4 . Then we set sin a = s 2 above, and also w = ssin# to obtain the 
Legendre normal form: 

/ g(z)dz = 2c(d)s 1 = (l + y RkSkl ) d9 

J-y Jo VI - sin 2 a sin 2 6 \ iTi 1 ~ Uk sm J 

This is a combination of one complete elliptic integral of the first kind and 4 
complete elliptic integrals of the third kind. The arguments n k — ~s 2 s~ k 2 of the 
integrals of the third kind are given in Table Q] 

Therefore, 

J g(z)dz = 2c(d) \ K{m) + ^ R k s k iIL(n k ; m)\ 

where K and II denote the complete elliptic integrals of the first and third kind, 
respectively. 

ASYMPTOTICS: s -»■ V2- 1, so m -> 0.029437251, a -> 0.172425997 rad 
~ 9°52'45.42". Also, si,s 2 -»• 1 and s 3 = 1 =-> (1 - \/2)/(l + v 7 ^). 

EXPERIMENTAL DATA: The hypergeometric functions were evaluated us- 
ing Romberg iteration. Coefficients and residues obtained symbolically and then 
numerically. Digits are not guaranteed to be all significative. 
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LI 


Err 1 Vd - 2 


Q 
o 


— n 9sm 34 


A 

■4: 


— n 31 Q97Q 


K 
o 






— 3480(14 


7 




8 


-0.360010 


9 


-0.363712 


10 


-0.366583 


10 2 


-0.387335 


10 3 


-0.389199 


10 4 


-0.389384 


10 5 


-0.389402 


10 6 


-0.389404 


10 7 


-0.389405 


10 8 


-0.389405 


10 9 


-0.389405 



Remark 4. Mark Rybowicz [15] provided the following alternative formula for 



Cd,i — — 



+ 



4d(u - 2) 



\fu{u — l)(u — d)ir 



K(v) 



u + 1 -n(-<^,t,) 



\fu{u — l)ir Au 
(2-d)(u + d) (d-u) 2 
y/u(d — u)ir 4du 



where 



2d , V2^m 
and u = 



d - 1 ~~ " 2 
His formula agrees with ours up to six decimal places. 



8 Appendix: Asymptotics for P n 

Let us recall briefly the large deviation principle for large radom matrices. A 
good reference is Guionnet [5] where we have taken most of the following. 

Let X be a real nxn symmetric matrix. Its entries are independent Gaussian 
random variables with mean and variances 2/n for a diagonal entry and 1/n 
for an off-diagonal one. Thus \f^X is in GOE(n). 

Let 'P(R) denotes the set of probability measures on R equipped with the 
weak topology. For any A e R" let (5(A) = ^ Y^i=i be the probability measure 
on R defined by 

6(X)(A) = -#{i : A t 6 A} . 
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We also need the function / : V(M) — ► [0, oo] defined by 

I (f*) = \J x 2 d(J l (x) -\ J \og{\x - y\)dn{x)dn{y) - | 
and the probability measure Q n on R™ with density 

® n = h II lAi-A.lexpf-^AfUA 

l<i<j<n \ i=l / 

with respect to the Lebesgue measure on K™ (Z™ is the normalisation constant). 
Then, according to Guionnet [8] Theorem 3.1, for every closed set A C V(M) 
one has 

Urn sup \ \ogQ n ({A e K" : 8(X) G .4}) < — inf 

To obtain an estimate on P n we take 

A = {n^V(R) : A»([0,oo[) = 1} 

so that, when A is the vector of eigenvalues of the matrix X, we have fj, E A 
if and only if X is semi-positive definite or (almost surely) if and only if X is 
positive definite. Since Q n is the joint law of the eigenvalues of X we get 

limsup — logProb{Jf £ S n : X is pd} < — inf 

where the probability is taken in ^J^GOE{n). Since the set of positive definite 
matrices is invariant under scaling it is not too difficult to see that 

P n = Prob {X e S n : X is pd} 

so that 

limsup — log(P n ) < — inf I(fi) = —a. 

It remains to explain why a is a positive number. The map I is > 0, each 
sub-level set {fi : I(fi) < M} is compact and, it achieves its minimum value at 
a unique probability measure on K described as Wigncr's semicircular law 

— sj 4 — x 2 dx 

with support at [—2,2]. Since this measure is clearly not in our set A we get 
a > and we are done. 
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